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Abstract:- alternating direction explicit and alternating direction implicit methods (ADE and ADI) were used 
to solve Schnakenberg model, we were found that alternating direction implicit method is much more accurate 
and faster than alternating direction explicit in this kind of models. 
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I. INTRODUCTION 

Reaction-diffusion (RD) systems arise frequently in the study of chemical and biological phenomena 
and are naturally modeled by parabolic partial differential equations (PDEs). The dynamics of RD systems has 
been the subject of intense research activity over the past decades. The reason is that RD system exhibit very 
rich dynamic behavior including periodic and quasi -periodic solutions [6]. 

Various orders are self-organized far from the chemical equilibrium. The theoretical procedures and 
notions to describe the dynamics of patterns formation have been developed for the last three decades [4]. 
Attempts have also been made to understand morphological orders in biology [5], Clarification of the 
mechanisms of the formation of orders and the relationship among them has been one of the fundamental 
problems in non-equilibrium statistical physics [3]. 

II. MATHEMATICAL MODEL 

A general class of nonlinear-diffusion system is in the form 



d ± Au + a ± u + b t v + f(u,v) + g ± (x) 

(1) 



dt 

dv 

— = d 2 Au + a 2 u + b 2 v - f(u, v) + g 2 (x) 

With homogenous Dirchlet or Neumann boundary condition on a bounded domain Q. , n<3, with 
locally Lipschitz continuous boundary. It is well known that reaction and diffusion of chemical or biochemical 
species can produce a variety of spatial patterns. This class of reaction diffusion systems includes some 
significant pattern formation equations arising from the modeling of kinetics of chemical or biochemical 
reactions and from the biological pattern formation theory. 

Schnackenberg model: 

ay = -k,b t = a 2 = b 2 = 0,/ = u 2 v, gi = a,g 2 = b, 
where k, a and b are positive constants. 

Then one obtains the following system of two nonlinearly coupled reaction-diffusion equations (the 
Schnakenberg model), 

— = d x hu -ku + u 2 v + a, (t, x) £ (0, <x>) x Q.\ 

dv 

— = d 2 Av-u 2 v + b, (t, x) £ (0, oo) x Q } 
at 

u(t,x) = v(t,x) = 0, t>0, xedQ. 
u(0.jO = UoGO, v(0,x) = i7 (x), xeQ ) 
Where d v d 2 , a, k and b are positive constants [8]. Various finite difference algorithms or schemes 
have been presented for the solution of hyperbolic -parabolic problem or its simpler derivatives, such as the 
classical diffusion equation. It is well-known that many of these schemes are partially unsatisfactory due to the 
formation of oscillations and numerical diffusion within the solution [1, 7]. 

Solution by the finite difference method, although more general, will involve stability and convergence 
problems, may require special handling of boundary conditions, and may require large computer storage and 
execution time. The problem of numerical dispersion for finite difference solutions is also difficult to overcome 
[2]. 
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DERIVATION OF ALTERNATING DIRECTION EXPLICIT (ADE) FOR Schnakenberg MODEL: 

The two dimensional Schnakenberg model is given by 

du (d 2 u d 2 u\ „ "\ 

— = dA —^ + —^ - ku + u 2 v + a, (t, x) £ (0, oo) x Q.\ 

dt \dx l dy l J " 

dv (d 2 v d 2 v\ 

Yt- d \^ 2+ W^ UV + h ' fc*)e(o.«)xn 



(2) 



We consider a square region 0<x<l, 0<y<l and u, v are known at all points within and on the 
boundary of the square region. We draw lines parallel to x, y, f-axis as x=\h, y=]k, and t=nz, ij=0, 1, 2,...,M 
and n=0, 1,2,..., N, where h=dx, k= 8y, z= St. 

The explicit finite difference approximation (ADE) to shnakenberg model in two-dimensions are given by: 
Uiin+i-Uiin (Ui-i in -2U iin + U i+lin \ fU u _ ln -2U iin + U u+ln \ , ,2 

+ a 

Viin+l-Viin (Vi- 1 , n -2V ijn +V i+1 , n \ (Vu- ln -2V ljn + V lj+ln \ . -2 

^v^ 1 = d2 ( w — ) + d2 ( w ) - ^ ;<^ + b 

Multiplying both equations by St and set h — k, then we have a rectangular region and replacing -^- = r x and 

— j- = r 2 ,then we get 

IWl - ^,;,n = 

ri{Ui- 1Jin - 2U iJin + Ui+lj, n + ty,;-l,n - 2£/;, ; -, n + £/y+i,„) - kSt U lJin + 

St (U lJin ) 2 V iJin + aSt 
and 

Vijfi+l - V lj,n = r 2(yi-l,j,n - 2V iJi7l + V i+lijin + Vij-xm ~ 2V lj,n + ^j+l,n) ~ St (£/;,;,„)%,„ + St b 

Then simplifying the system to obtain 

Ui,j, n +i = (1 - 4fi - kSt)U iJin + r x {Ui- t j in + U i+iihn + U u _ lin +U u+lin ) + St (£/;,,>)%,„ + aSt 
Vi,j, n+ i = (1 " 4r 2 )V y , n + r 2 (V £ _ v , n + K i+1 , y , n + Vy.^ + V lJ+lin ) - St (t/y, n )Vy, n + St b 
This is the alternating direction explicit formula for the Schnakenberg model 

III. DERIVATION OF ALTERNATING DIRECTION IMPLICIT (ADI) FOR 

GLYCOLYSIS MODEL 

In the ADI approach, the finite difference equations are written in terms of quantities at two x levels. 
However, two different finite difference approximations are used alternately, one to advance the calculations 
from the plane n to a plane n+1, and the second to advance the calculations from (n+7)-plane to the (n+2)- 
plane by replacing — -y and — j by implicit finite difference approximation [5]. we get 

si ~ dl \ h 2 ) + dl \ k 2 ) ~ kUi ^ n 

+ {y iijin ) 2 v iijin + a 

Vtj * +l s ; Vtj * = d 2 ( v i-u»«- 2v u»« +v w»^ + d 2 ^-**-yiw) - {u iJin )\ jin + b 

we set h — k, then we have a rectangular region and multiplying the equations by St and replacing -^- = r x 

and — — = r 2 , then we get 

Uij, n+ i - £/y,„ = r 1 (f/i_ V) „ +1 - 2U lJ>n+1 + f/ i+ i,;,„ + i + Ui,j-i, n ~ 2U lJiTl + U iJ+ltn ) - kSt U lJin 

+ St(U i j in ) 2 V i j n + Sta 
Vij, n+ i - V tJ , n = r 2 (y t - tJ , n+1 - 2V t j, n+1 + V i+tJ , n+t + Vu-^ - 2Vy, n + V iJ+t , n ) - St (l/y,„) 2 V y , n + St b 

-r2Vi-ij,n + i + (1 - 2r 2 )V iJin+1 - r 2 V i+1J +1 , n+1 

= r 2 Vi,j-i,n + (1 - 2r 2 )V iJin + r 2 V lJ+lin - St (U tJfl ) 2 V iJjn + Stb 
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-riUi-ij,n+i + (1 - 2r t )U iJin+ i - rt/ f+1<;>+1 = (1 - 2r x )U iijin + r t U Uj+tin + Sta + 8t (U iijin ) 2 V iiJin 
or 
-ri£/i-i, M+ i + (1 - 2r 1 )f/ iJ , n+1 - r x U i+UA+1 

= (1 - 2r 1 )£/y,„ + nUij+^n + St (f/y,„)Vy, n + 5t (a - fct/y,„) 
and 
-r 2 ^-ij,n + i + (1 - 2r 2 )lA ;>+1 - r 2 V f+1J+1<n+1 

= r 2 Vij_ lin + (1 - 2r 2 )Vy )n + r 2 Vy +1 , n - St (£/y,„)Vy, n + 5t fot 
Advance the solution from the (n+l) th plane to (n+2) th plane by replacing — \ and — \ with explicit finite 
difference approximation at (n+l) th plane then — | and — \ by an implicit finite approximation at the (n+2) th 
plane 

"lJ.n+2 ~ UlJ. n + l , ( Ul-lJ.n + l ~ 2U tJ , n+1 + ^ +1 ,;,n +1 \ , . ( UjJ-l.n+2 ~ ^i,j,n + 2 + UjJ + l.n +2 \ 

si = dl I # J + dl I P J 

-fcl/y,„ + (X;,n)\;,n + a 

and 

^-^ - * C--"- 2 ";?" -- ) + d. r"--- 2 ^r"""-" ) - (»„,„) v„ + » 

Multiplying the equations by St and replacing r x — -^- and r 2 = -^- when h — k, then we have 

^,;,n+2 - U iJin+1 = r^Ut-y^i ~ 2[/y,„ +1 + U i+1J „ +1 ) + n(U \j _ 1>n+2 - 2(/y >+2 + £/y + l,„ +2 ) 

- fc5t f/y <n + 5t (Uij ,n)\),n +8ta 
and 

»W+2 - V iJin+1 = ^(Vi-y.n+i - 2V t j, n+1 + V l+1J , n+1 ) + r 2 (7y_!, n+2 - 27y >+2 + 7y +1 ,„ +2 ) 

- 5t (f/y,„) Vy, n + 5tfo 

and this implies that 

-rif/y-i.n+2 + (1 + 2r 1 ){/y >+2 - nUij+^+2 

= nUi-ij,n + (1 - 2r 1 )f/y j+ln + r 1 f/ i+lj + n+1 + 5t (a - fc£/y>) + (tfij>)Vi,;,n 
and 
-r2^,;-l,n + 2 + (1 - 2r 2 )V iJin+2 - r 2 7y +1 ,„ +2 

= ^-ij ,„+i + (1 - 2r 2 )7y,„ +1 + r 2 7 i+y , n+1 - 5t (f/y, n )Vy, n + St b 
Then we get two systems 

-nU t ,J -l, n +2 + (1 + 2r 1 )f/y,„ +2 - rjt/y+j^+j 

= r l^-lj,n + (1 - 2r 1 )f/y,+i„ + r 1 f/ i+lj + n+1 + St (a - fcf/y,,) + (£/y>)Vy,„ 
-^2^;-i,n+2 + (l-2r 2 )7y, n+2 -r 2 V'y +1< „' +2 

= r 2 ^_ 1J , n+1 + (1 - 2r 2 )Fy, n+1 + r 2 V mj>+1 - St {U iJin )\j, n + St b 
and 
-riU t - tJin+1 + (1 - 2r 1 )f/y,„ +1 - ri£/ f+1<; - n+1 

= (1 - 2r 1 )f/y > + r l U lij+lin + St (£/y,„) 2 Vy, n + ft (a - fcl/y,„) 
-^V i -y, n+ i + (l-2r 2 )Vy, n+1 -r 2 7 mj+1 ,l +1 

= r 2 Fy_i,„ + (1 - 2r 2 )Vy <n + r 2 7y+i, n - ft (£/y,„)Vy,„ + ft ft 

The last two systems represent Alternating Direction implicit method under the conditions 

u l,q,n+l = u l,q + l,n+l =° 
u M,q,n + \ = U M ,q + l,n + \ = °- 

and 

V l,q,n+\ = V l,q+l,n+l = " 
V M , q ,n + l =V M ^ J+Ul+1 =0. 

The tridiagonal matrices for the system in the level n advanced to the level n+l, for both u and v can be 
formulated as follows AU=B. 
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(1+2*1) -n 

-n (1 + 2,)) 



-r, (l + 2rj) - r x 

-r\ (l + 2rj) 



"4, 9 ,n+l 

"5.</+l 



u ii-3, 9 +1 

"n-2,?+l 

_"M-l, ? ,n+l 



"Z + r 2 ( u 2,q+l,n + u 2,q-l,n ) + (1 ~ 2r 2 - z)« 2 , 9 ,n + Z« ^A,n v 2,q,n 
<**■ + r 2 ("3,?+l,n + M 3,g-l,n ) + (1 - 2r 2 - z)«3 j?! „ + z« 3, ? ,nv 3)?>n 
az + r 2 (»4, 9 +l,n + u A,q-l,n ) + (1 - 2r 2 - z)u 4qn + zu A,q,nv Aqn 



"z+r 2 (u M _ lq+ln +u M _ lq _ ln ) + (l-2r 2 -z)u M _ lqn + zu M-l,q,nv M _ lqn 



(1 + 2))) -r { . . 

-/) (l + 2)j) -/j 

-/) (1 + 2^) -rj 

-n (l + 2n) -r. 



-,i (1 + 2,]) -/j 
-IJ (l + 2jj) 



v 3,<7,n+l 
l '5, 9 ,n+l 



"M-3,?,n+l 
l 'M-2,g,rc+l 
v iW-l, g ,n+l _ 



''2( l '2, 9 +l,« + v 2,«-l ; n) + (! - 2, 2 - - K ) v 2,q,n - » Ul^l.q.n 
r 2( v 3,q+ln + v 3, ? -l,« ) + C 1 " 2, 2 " ^ v 3, ? ,n " ZM 3,?.» v 3, 9 ,« 
^( v 4,9+l,« + v 4, 9 -l,n) + ( ! " 2, 2 " Z K )%q,n " Z" 2 4,?.« v 4, ? ,n 



.) + (1-2^-^4 



And the tridiagonal 
given by AU=B. 



for the system in level «+l advanced to level n+2 for both u and v are 
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u p,2 


,n+2 


l 'p,3,n+2 


"pA 


,n+2 


U P,5 


,«+2 


"p,M 


-3,n+2 


"p,M 


-2,77+2 


"p.M 


-l,n+2 _ 



2 '2> -'2 ° ° 

r, (l + 2r 2 ) -r 2 

-q (1 + 2/!) -i) 

-r 2 (l + 2,' 2 ) 



-i, (l + 2r 2 ) 

PZ + l(«p+l,2,n+l +"/?-l,2,n+l) + ( 1 ^ 2 1 ~ z )" p,2,n+\ ~ zu p,2,n + ^ P, 2 . nV p,2,n 
fz + n ( B p+1,3,»+1 + "/>-l,3,n+l ) + (' - 2r l - z)»p,3,n+l " z "/j,3,h + z " P,3,'» v p,3,n 
P z + l("/;+l,4,«+l +'V-l,4,«+l> + < 1 -2'- 1 -z)« pA „+i -z«p,4,« + z " pA»»p,4,» 
Z^ + 'l^+U.n+l +M i)-l,5,«+l ) + (1_2 ''l _z),< p,5,n+l ~ z V5,« + zu />,5,n>,5,n 



l,n+l + "p-l,M-l,«+l) + (1 - 2q - :)« p? M-l,n+l " zu pM-\,2,n + z " 

Also for AV=B the tridiagonal is in the form 



(l + 2m 2 ) -m 2 0. 

-m 2 (l + 2m 2 ) -m 2 

-m 2 (l + 2ra 2 ) -m 2 

-m 2 (l + 2m 2 ) -'» 2 

-mj 



(1 + 2m 2 ) - m 2 



p.l.n-2 
p,\n+2 
pA-n+2 
p,5,n+2 



v p,M~3,n+2 
v p,M-2,n+2 
v p,M-\,n+2 



m l( v p+l,2,n+l +v p-l,2,n+0 + ( - l ' 2m l^ zK '> v p,2,n+'i +zu p,2,n~ zu P^ nV p,2,n 
m l0>+l,3,n+l +v p-l,3,n+l) + ( 1 - 2m l ~ zK >p2i,n+\ + zu p,3,n ~ zu pXnVp^n 
'"1 (>+l,4,n+l + v p-l,4,«+l) + (! " 2m l " zK ~> v pA,n+l + -'"p,4,n " *" pA,n v p A,n 



'1 ( v p+l,M-l,n+l + v p-l,jtf -l,n+l ) + 



> v p,M-ln+l + zu p,M-ln ~ zu 1 



IV. Application (numerical example) 

We solved the following example numerically to illustrate the efficiency of the presented methods, 
suppose we have the system 
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-2u + u 2 v + a, 
-u 2 v + u + b, 



dt 



(d 2 u d 2 u\ 
(d 2 v d 2 v\ 



We the initial conditions 

u(x, y, 0), v(x,y, 0)) = (exp(-x - y) , expU + y)) 
u(0,y,t),K0,y,O) = (exp(-t/2-y),exp(t/2+y)) 
u(x, 0, t), v{x, 0, t)) = (exp(-t/2 - x) , exp(t/2 + x)) 

(u, vf) — (exp i x — y) , exp [- + x + y)) is the exact solution of the problem. 

Where a=b=0.1. 

Then the results in more details are shown in following table and figure: 




X 


ADE 


ADI 


Exact 


0.1 


2.69265 


2.6912344 


2.6912344 


0.15 


2.423622 


2.27384012 


2.20573840 


0.2 


2.181472 


2.01678803 


1.8971263 


0.25 


1.963516 


1.8279602 


1.6823796 


0.3 


1.767336 


1.64066 


1.52614066 


0.35 


1.590757 


1.4899208 


1.401599208 


0.4 


1.431821 


1.31680 


1.254730168 


0.45 


1.288764 


1.18920286 


1.14620286 


0.5 


1.160001 


1.0435533 


1.011435533 


0.55 


1.044102 


0.939289 


0.879392894 


0.6 


0.939783 


0.8454428 


0.78454428 


0.65 


0.845887 


0.7609726 


0.697609726 


0.7 


0.761373 


0.6849420 


0.635849420 


0.75 


0.685302 


0.6165078 


0.568960165 


0.8 


0.616832 


0.5549110 


0.52549110 


0.85 


0.5552032 


0.4994686 


0.44994686 


0.9 


0.4997315 


0.4495655 


0.40495655 


0.95 


0.4498021 


0.4046483 


0.374046483 


1 


0.4048613 


0.403642189 


0.363642189 



We saw that alternating di 
for solving Schnakenberg model. 



impli 



CONCLUSION 

it is more accurate than alternating direction explicit method 
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